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Abstract The famous three-body problem is investigated by means of a numerical 
approach with negligible numerical noises in a long enough time interval, namely the 
Clean Numerical Simulation (CNS). From physical viewpoints, position of any bodies 
contains inherent micro-level uncertainty. The evaluations of such kind of inherent 
\*J micro-level uncertainty are accurately simulated by means of the CNS. Our reliable, 

very accurate CNS results indicate that the inherent micro-level uncertainty of position 
of a star/planet might transfer into macroscopic randomness. Thus, the inherent 
micro-level uncertainty of a body might be an origin of macroscopic randomness of 
the universe. In addition, from physical viewpoints, orbits of some three-body systems 
at large time are inherently random, and thus it has no physical meanings to talk 
about the accurate long-term prediction of the chaotic orbits. Note that such kind of 
uncertainty and randomness has nothing to do with the ability of human being. All of 
these might enrich our knowledge and deepen our understandings about not only the 
three-body problem but also chaos. 

Key Words Three-body problem, chaos, multiple precision, Taylor expansion, micro- 
level uncertainty 
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1 Introduction 

It is well-known that the microscopic phenomena are inherent random, while many 
macroscopic phenomena such as moving stars and planets in the universe looks ran- 
dom as well. What is the origin of the macroscopic randomness? Are there any 
relationships between the microscopic uncertainty and the macroscopic randomness? 
In this article, using chaotic motion of the famous three-body problem as an example, 
we illustrate that the micro-level uncertainty of position of stars/planets might be one 
origin of the macroscopic randomness of the universe. 

It is a common knowledge that some "deterministic" dynamic systems have chaotic 
property: their numerical simulations have sensitive dependence on initial conditions 
(SDIC), i.e. the so-called butterfly-effect, so that long-term accurate prediction is 
impossible (8 10 . Since truncation and round-off errors are unavailable for all numer- 



ical simulation techniques, nearly all numerical results given by traditional methods 



based on double precision are not "clean" : they are something mixd with the so-called 
"numerical noises". Due to the SDIC, truncation and round-off errors enlarge expo- 
nentially so that it is very hard to gain reliable chaotic results in a long time interval. 
As pointed out by Lorenz n0\ in 2006, different traditional numerical schemes (based 
on 16 or 32 digit precision) may lead to not only the uncertainty in prediction but 
also fundamentally different regimes of solution. 



In order to gain reliable chaotic results in a long enough time interval, Liao [5] 
developed a numerical technique with negligible numerical noises, called the "Clean 
Numerical Simulation" (CNS). Using the computer algebra system Mathematica with 
the 400th-order Taylor expansion and data in 480-digit precision, Liao [5] obtained, 
for the first time, the reliable numerical results of chaotic solution of Lorenz equation 
in a long time interval < t < 1000 Lorenz time unit (LTU). Liao's "clean" chaotic 
solution of Lorenz equation was confirmed by Wang et al. |17], who employed the 
parallel computation and the multiple precision (MP) library to gain reliable chaotic 
solution up to 2500 LTU by means of the CNS approach with lOOOth-order Taylor 
expansion and data in 2100-digit precision, and their result agrees well with Liao's 
one (5J in < t < 1000 LTU. This confirms the validity of the CNS approach. 

It was found by Liao [5] that, to gain a reliable "clean" chaotic solution of Lorenz 
equation in the interval < t < T c , the initial conditions must be at least in the 
accuracy of 10 _2Tc//5 . For example, in the case of T c = 1000 LTU, the initial condition 
must be in 400-digit precision at least. It should be emphasized that, the 400-digit 
precision, which is "mathematically" necessary for the initial condition and all data 
at each time-step, is so high that even the statistical fluctuation of velocity and tem- 
perature becomes a very important physical factor and therefore cannot be neglected. 
However, from the physical viewpoints, Lorenz equation (as a macroscopical model 
for climate prediction on Earth) completely neglects the influence of the statistic fluc- 
tuation of velocity and temperature about the climate. Therefore, as pointed out by 
Liao k)J , this leads to the so-called "precision paradox of chaos" . 

How to avoid such kind of paradox? Traditionally, it is believed that Lorenz equa- 
tion is a "deterministic" system, say, its initial condition and all physical parameters 
are completely certain, i.e. "absolutely accurate" . However, such kind of "absolutely 
accurate" variables only exist in mathematics, which have no physical meanings in 
practice, as mentioned below. For example, velocity and temperature of fluid are 
concepts defined by statistics. It is well-known that any statistical variables con- 
tain statistic fluctuation. So, strictly speaking, velocity and temperature of fluid are 
not "absolutely accurate" in physics. Using Lorenz equation as an example, Liao [6] 
pointed out that its initial conditions have fluctuations in the micro-level 10 -30 so that 
Lorenz equation is not deterministic, from the physical viewpoint. Although 10~ 30 is 
much smaller than truncation and round-off errors of traditional numerical methods 
based on double precision, it is much larger than 10~ 400 that can be used in the CNS 
approach. Thus, by means of the CNS, Liao (6) accurately simulated the evaluation of 
the micro-level uncertainty of initial condition of Lorenz equation, and found that the 
micro-level uncertainty transfers into the observable randomness. Therefore, chaos 
might be a bridge between the micro-level uncertainty and the macroscopic random- 



ness, as pointed out by Liao |6]. Currently, Liao [7] employed the CNS to the chaotic 
Hamiltonian Henon-Heiles system for motion of stars orbiting in a plane about the 
galactic center, and confirmed that, due to the SDIC, the inherent micro-level uncer- 
tainty of position of stars indeed evaluates into the macroscopic randomness. 

However, Lorenz equation is a greatly simplified model of Navier-Stokes equation 
for flows of fluid. Besides, unlike Hamiltonian Henon-Heiles system for motion of 
stars orbiting in a plane, orbits of stars are three dimensional in practice. Thus, in 
order to further confirm the above conclusion, it is necessary to investigate some more 



accurate physical models, such as the famous three-body problem 111 3, 16 governed 



by the Newtonian gravitation law. In fact, non-periodic results were first found by 
Poincare [II] for three-body problem. In this paper, using the three-body problem as 
a better physical model, we employ the CNS to confirm the conclusion: the inherent 
micro-level uncertainty of position of a star /planet might transfer into the observable, 
macroscopic randomness of its orbit so that the inherent micro-level uncertainty of a 
star /planet might be an origin of the macroscopic randomness of the universe. 



2 Approach of Clean Numerical Simulation 

Let us consider the famous three-body problem, say, the motion of three celestial 
bodies under their mutual gravitational attraction. Let Xi, Xi, £3 denote the three or- 
thogonal axises. The position vector of the i body is expressed by r^ = (xij, X2,i, £3,1). 
Let T and L denote the characteristic time and length scales, and mi the mass of 
the zth body, respectively. Using Newtonian gravitation law, the motion of the three 
bodies are governed by the corresponding non-dimensional equations 
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denotes the ratio of the mass. 

In the frame of the CNS, we use the M-order Taylor expansion 
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to accurately calculate the orbits of the three bodies, where the coefficient a£f is only 
dependent upon the time t . Note that the position Xk,i{t) and velocity &&,»(£) at 



t = to are known, i.e. 

«o' 1 = Xk,i(to), «/ = Xk,i(to). (5) 

The recursion formula of a%£ for m > 2 is derived from Q, as described below. 

Write 1/Rij in the Taylor expansion 

M 

/« = #T « E ^ (* - *>)* ( 6 ) 

with the symmetry property j3%? = (3^, where (3^ is determined later. Substituting 
Q and (J6| into ([!]) and comparing the like-power of (t — t ), we have the recursion 
formula 
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Thus, the positions and velocities of the three bodies at the next time-step to + At 
read 

M 

x k>i (to + At) « ^^(At)™, (8) 

m=0 

M 

x k)i (t + At) « ^m^MAt)™- 1 . (9) 

m=l 

Write 

M M 

^ = <• « E t# (* - *o) m > & = E <#(* - t o) m ' ( 10 ) 
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with the symmetry property 7,^ = 7^ and cr^f = o{^. Substituting (J2J), Q and (J6J) 
into the above definitions and comparing the like-power of t — t , we have 
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with 
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tti = E E (^ ~ ^n) {°^~n ~ °^~r) > i?3,m>l, (13) 

fe=l n=0 

and the symmetry //^ = /i^\ Using the definition ([6]), we have 



Substituting (10) into the above equation and comparing the like-power of (t — to), 
we have 
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3 a^-n = 0, m > 1, 



n=0 

which gives the recursion formula 
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In addition, it is straightforward that 

/3o J = i J# = -Ry, c^=(#") 2 , 7o J = (^) 3 ' att = t - (15) 



i£ 



*../ 



It is a common knowledge that numerical methods always contain truncation and 
round-off errors. To decrease the round-off error, we express the positions, velocities, 
physical parameters and all related data in iV-digit precision, where N is a large 
enough positive integer. Obviously, the larger the value of N, the smaller the round- 
off error. Besides, the higher the order M of Taylor expansion pL the smaller the 
truncation error. Therefore, if the order M of Taylor expansion (El) is high enough 
and all data are expressed in accuracy of long enough digits, truncation and round-off 
errors of the above-mentioned CNS approach (with reasonable time step At) can be 
so small that numerical noises are negligible in a given (long enough) time interval, 
say, we can gain "clean", reliable chaotic numerical results with certainty in a given 
time interval. In this way, the orbits of the three bodies can be calculated, accurately 
and correctly, step by step. 

In this article, the computer algebra system Mathematica is employed. By means 
of the Mathematica, it is rather convenient to express all datas in 300-digit precision, 
i.e. N = 300. In this way, the round-off error is so small that it is almost negligible in 
the given time interval (0 < t < 1000). And the accuracy of the CNS results increases 
as the order M of Taylor expansion Q enlarges, as shown in the next section. 



3 A special example 

Without loss of generality, let us consider the motion of three bodies with the initial 
positions 

n = (5, 0, -1), r 2 = (0, 0, 0), r 3 = -(r x + r 2 ), (16) 

and the initial velocities 

ri = (0, -1, 0), r 2 = (1, 1, 0), r 3 = -fa + r 2 ), (17) 

where 6 is a constant. Note that S is the only one unknown parameter in the initial 
condition. For simplicity, we first only consider the three different cases: 5 = 0, 



5 = +10~ 60 and 5 = — 1CT 60 . Note that the initial velocities are the same in the 
three cases. Mathematically, the three initial positions have the tiny difference in the 
level of 10~ 60 , which however leads to huge difference of orbits of the three bodies at 
t = 1000, as shown below. For the sake of simplicity, let us consider the case of the 
three bodies with equal masses, i.e. pj — 1 (J — 1,2,3). We are interested in the 
orbits of the three bodies in the time interval < t < 1000. 

Note that the initial conditions satisfy 

3 3 

3=1 3=1 

Thus, due to the momentum conversation, we have 

3 3 

£>i(*) = 5>(*) = 0, *>0 (18) 

3=1 3=1 

in general. 

All data are expressed in 300-digit precision, i.e. N = 300. Thus, the round-off 
error is almost negligible. In addition, the higher the order M of Taylor expansion 
(El), the smaller the truncation error, i.e. the more accurate the results at t = 1000. 
Assume that, at t = 1000, we have the result xi : i = 1.8151012345 by means of the 
Mi-order Taylor expansion and the result xi 5 i = 1.8151047535 by means of the M2- 
order Taylor expansion, respectively, where M2 > M\. Then, the result xi,i by means 
of the lower-order (Mi) Taylor expansion is said to be in the accuracy of 5 significance 
digit, expressed by n s = 5. For more details about the CNS, please refer to Liao [7]. 

When 5 = 0, the corresponding three-body problem has chaotic orbits with the 
Lyapunov exponent A = 0.1681, as pointed by Sprott [13] (see Figure 6.15 on page 
137). It is well-known that a chaotic dynamic system has the sensitivity dependence 
on initial condition ( SDIC). Thus, in order to gain reliable numerical results of the 
chaotic orbits in such a long interval < t < 1000, we employ the CNS approach 
using the high-enough order M of Taylor expansion with all data expressed in 300-digit 
precision. 

It is found that, when 5 = 0, the CNS results at t = 1000 by means of At = 10~ 2 , 
N = 300 and M = 8, 16,24,30,40 and 50 agree each other in the accuracy of 11, 
23, 38, 48, 64 and 81 significance digits, respectively. Approximately, n s , the number 
of significance digits of the positions at t = 1000, is linearly proportional to M (the 
order of Taylor expansion), say, n s ~ 1.6762M — 2.7662, as shown in Fig. [I] For 
example, the CNS approach using the 50th-order Taylor expansion and data in 300- 
digit precision (with At = 10~ 2 ) provides us the position of Body 1 at t = 1000 in 




Figure 1: Accuracy (expressed by n s , number of significance digits) of the CNS 
results at t = 1000 in the case of 5 = by means of N = 300, the different time-step 
At and the different order (M) of Taylor expansion. Square: At = 10~ 2 ; Circle: 
At = 10~ 3 . Solid line: n, « 1.6762M - 2.7662; dashed line: n, « 2.7182M - 4.2546. 



the accuracy of 81 significance digit: 



xi,i = +1.8151047535629516172165940088454400645690 
03032055743237590103285240436181354986834, 

x 2) i = -1.4406351440582861611338350554067489001231 
29002848537219768908176322703247482409938, 

x 3 ,i = +1.9870078875786298810976776414985778979168 
46299746397051757074.117730821734512844759. 



(19) 
(20) 
(21) 



Using the smaller time step A = 10~ 3 and data in 300-digit precision (i.e. N = 300), 
the CNS results given by the 8, 16, 24 and 30th-order Taylor expansion agree in the 
accuracy of 18, 38, 62 and 77 significance digits, respectively. Approximately, n s , the 
number of significance digits of the positions at t = 1000, is linearly proportional to 
M (the order of Taylor expansion), say, n s « 2.7182M — 4.2546, as shown in Fig. [I] 
It should be emphasized that the CNS results by At = 10~ 3 and M = 30 agree 
(at least) in the 77 significance digits with those by At = 10~ 2 and M = 50 in the 



whole time interval < t < 1000. In addition, the momentum conservation (18) is 
satisfied in the level of 10 -295 . All of these confirm the correction and reliableness 
of our CNS result^] Thus, although the considered three-body problem has chaotic 
orbits, our numerical results given by the CNS using the 50th-order Taylor expansion 
and accurate data in 300-digit precision (with At = 10~ 2 ) are reliable in the accuracy 
of 77 significance digits in the whole interval < t < 1000. 



*Liao u\ proved a convergence-theorem and explained the validity and reasonableness of the CNS 
by using the mapping x n +\ = mod(2a;„, 1). 





Figure 2: x — y and x — z of Body 1 (0 < t < 1000) in the case of 5 = 0. 




Figure 3: x — y and x — z of Body 2 (0 < £ < 1000) in the case of 6 = 0. 
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Figure 4: x — y and x — z of Body 3 (0 < t < 1000) in the case of 5 = 0. 
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Figure 5: Accuracy (expressed by n s , number of significance digits) of the result at 
t = 1000 in the case of 5 = 10 -60 by means of N = 300, the different time-step At and 
the different order of Taylor expansion (M). Square: At = 10~ 2 ; Circle: At = 10 -3 . 
Solid line: n, « 0.2885M - 3.4684; dashed line: n, « 1.5386M - 3.7472. 



The orbits of the three bodies in the case of 8 = are as shown in Figs. [2] to 
|1J The orbits of Body 1 and Body 3 are chaotic. This agrees well with Sprott's 
conclusion [13] (see Figure 6.15 on page 137). However, it is interesting that Body 2 



oscillates along a line on the plane z — 0. So, since ^, =1 ij = Ylj=i r j = due to the 
momentum conservation, the chaotic orbits of Body 1 and Body 3 must be symmetric 
about the regular orbit of Body 2. Thus, although the orbits of Body 1 and Body 3 
are disorderly, the three bodies as a system have an elegant structure with symmetry. 

The initial conditions when 8 = 10~ 60 have a tiny difference 

dr 1 = 10" 60 (1,0,0) 

from those when 5 = 0. Thus, it is reasonable to assume that the corresponding 
dynamic system is chaotic, too. Similarly, the corresponding orbits of the three bodies 
can be accurately simulated by means of the CNS. It is found that, when 5 = 10 -60 , 
the CNS results at t = 1000 by means of At = 10" 2 , N = 300 and M = 16, 24, 
30, 40, 50, 60 , 70, 80, 100 agree each other in the accuracy of 1, 3, 5, 9, 11, 14, 
17, 19 and 27 significance digits, respectively. Approximately, n s (the number of 
significance digits) is linearly proportional to M (the order of Taylor expansion), say, 
n s w 0.2885M — 3.4684, as shown in Fig. [5} According to this formula, in order to 
have the CNS results (at t = 1000) in the precision of 81 significance digits by means 
of At = 10~ 2 , the 300th-order of Taylor expansion, i.e. M = 300, must be used. But, 
this needs much more CPU time. To confirm the correction of these CNS results, we 
further use the smaller time step At = 10 -3 . It is found that, when 5 = 10~ 60 , the 
results at t = 1000 by means of the CNS using At = 10" 3 , iV = 300 and M = 8, 16, 
24, 30, 40, 50 agree well in the precision of 8, 21, 33, 43, 59, 72 significance digits, 
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Figure 6: x — y and x — z of Body 1 (0 < t < 1000) in the case of S = 10 
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Figure 7: x — y and x — z of Body 2 (0 < t < 1000) in the case of S = 10 
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Figure 8: x — y and x — z of Body 3 (0 < t < 1000) in the case of S = 10 
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Figure 9: Orbit of Body 1 (0 < t < 1000). Left: 5 = 0; Right: 5 = 10 
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Figure 10: Orbit of Body 2 (0 < t < 1000). Left: 5 = 0; Right: 5 = 10 
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Figure 11: Orbit of Body 3 (0 < t < 1000). Left: 5 = 0; Right: 5 = 10 
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respectively. Approximately, n s , the number of significance digits of the corresponding 
results at t — 1000, is linearly proportional to M (the order of Taylor expansion), say, 
n s m 1.5386M — 3.7472, as shown in Fig. [5] For example, the position of Body 1 at 
t = 1000 given by the 50th-order Taylor expansion and data in 300-digit precision 
with At = 10~ 3 reads 



Xll = +10.5718991771626848605399651180233387355185816 

19035776331820966524368101464, (22) 

x 2>1 = -33.3956860196582147031781512361768602407559680 

29199276186710038142873929480, (23) 

x 3>1 = +20.2845527396821922952136869217938844104404153 

94341528671055848805097262214, (24) 

which are in the precision of 72 significance digits. Note that the positions of the 
three bodies at t = 1000 given by At = 10~ 2 and M = 100 agree well (in precision 
of 27 significance digits) with those by At = 10~ 3 and M = 50. In addition, the 



momentum conservation (18) is satisfied in the level of 10 293 . Thus, our CNS results 
in the case of 5 = 10~ 60 are reliable in the interval < t < 1000 as well. 

The orbits of the three bodies in the case of 5 = 10 -60 are as shown in Figs. [6J 
to |8j It is found that, in the time interval < t < 800, the orbits of the three 
bodies are not obviously different from those in the case of 5 = 0, say, Body 2 
oscillates along the same line on z = 0, Body 1 and Body 3 are chaotic with the 
same symmetry about the regular orbit of Body 2. However, the obvious difference 
of orbits appears when t > 810: Body 2 departs from the oscillations along the line 
on z = and escapes (together with Body 3) along a complicated three-dimensional 
orbit. Besides, Body 1 and Body 3 escape in the opposite direction without any 
symmetry. As shown in Figs. [9] to 11, Body 2 and Body 3 go far and far away from 



Body 1 and thus become a binary-body system. Thus, it is very interesting that, the 
tiny difference dv\ = 10 _60 (1,0, 0) of the initial conditions finally disrupts not only 
the elegant symmetry of the orbits but also even the three-body system itself! 

Similarly, in the case of 5 = — 10 -60 , we gain the reliable orbits of the three bodies 
by means of the CNS with At = 10" 3 , N = 300 and M = 20 (i.e. the 20th-order 



Taylor expansion). As shown in Figs. 12 to 14, the tiny difference in the initial 



position disrupts not only the elegant symmetry of the orbits but also the three- 
body system itself as well: when t > 810, Body 2 departs from its oscillation along 
the line on z = and escapes (but together with Body 1) in a complicated three- 
dimensional orbit, while Body 1 and Body 3 escape in the opposite direction without 
any symmetry. Note that, in the case of 5 = — 10 -60 , Body 1 and Body 2 go together 
far and far away from Body 3 to become a binary system. However, in the case of 
5 = +10~ 60 , Body 2 and Body 3 escape together to become a binary system! This is 
very interesting. Thus, the orbits of the three-body system when 5 = 0, 5 = +10 -60 
and 5 = — 10~ 60 are completely different. 

From the mathematical viewpoints, the above results are not surprising at all and 
thus there exist nothing new: since the three-body system is chaotic (as pointed out by 
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Figure 12: Orbit of Body 1 (0 < t < 1000). Left: 6 = +10" 60 ; Right: 5 = -10 
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Figure 13: Orbit of Body 2 (0 < t < 1000). Left: 5 = +10" 60 ; Right: 5 = -10 
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Figure 14: Orbit of Body 3 (0 < t < 1000). Left: 5 = +10^ 60 ; Right: 5 = -10 
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Sprott [13]), the results are certainly very sensitive to the initial conditions. However, 
the difference of the three initial positions is so small that they can be regarded as the 
same in physics! In other words, from physical viewpoints, such a small difference 
in space has no physical meanings at all, and thus the three initial conditions (when 
5 = and 5 = ±1CT 60 , respectively) are the same in physics. This is mainly because 
position of any a body inherently contains the micro-level uncertainty so that the 
three-body system is not deterministic, as explained below. 

It is well known that the microscopic phenomenon are essentially uncertain/random. 
Let us first consider some typical length scales of microscopic phenomenon which are 
widely used in modern physics. For example, Bohr radius 

h 2 

5.2917720859(36) x lO -11 (m) 



m P e 2 



is the approximate size of a hydrogen atom, where h is a reduced Planck's constant, 
m e is the electron mass, and e is the elementary charge, respectively. Besides, the 
so-called Planck length 



l p = J^f ~ 1.616252(81) x 1(T 35 (m) (25) 

is the length scale at which quantum mechanics, gravity and relativity all interact very 
strongly, where c is the speed of light in a vacuum and G is the gravitational constant. 



According to the string theory 12 , the Planck length is the order of magnitude of 
oscillating strings that form elementary particles, and shorter length do not make 
physical senses. Especially, in some forms of quantum gravity, it becomes impossible 
to determine the difference between two locations less than one Planck length apart. 
Therefore, in the level of the Planck length, position of a body is inherently uncertain. 
This kind of microscopic physical uncertainty is inherent and has nothing to do with 
Heisenberg uncertainty principle (4J and the ability of human being as well. 

In addition, according to de Broglie (2), any a body has the so-called wave-particle 
duality. The de Broglie's wave of a body has non-zero amplitude. Thus, position 
of a body is uncertain: it could be almost anywhere along de Broglie's wave packet. 
Thus, according to the de Broglie's wave-particle duality, position of a star/planet is 
inherent uncertain, too. Therefore, it is reasonable to assume that, from the physical 
viewpoint, the micro-level inherent fluctuation of position of a body shorter than the 
Planck length l p is essentially uncertain and/or random. 

To make the Planck length l p w 1.62 x 1CT 35 (m) dimensionless, we use the dimeter 
of Milky Way Galaxy as the characteristic length, say, du w 10 5 light year «9x 
10 20 meter. Obviously, l v jdu ~ 1.8 x 10~ 56 is a rather small dimensionless number. 
Thus, as mentioned above, two (dimensionless) positions shorter than 10~ 56 do not 
make physical senses in many cases. So, it is reasonable to assume that the inherent 
uncertainty of the dimensionless position of a star /planet is in the micro-level 10 -60 . 
Therefore, the tiny difference dx\ = ±10 _60 (1, 0, 0) of the initial conditions is in the 
micro-level: the difference is so small that all of these initial conditions can be regarded 
as the same in physics! 
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Mathematically, 10~ 60 is a tiny number, which is much smaller than truncation 
and round-off errors of traditional numerical approaches based on data of 16-digit 
precision. So, it is impossible to investigate the influence and evaluation of this inher- 
ent micro-level uncertainty of initial conditions by means of the traditional numerical 
approaches. However, the micro-level uncertainty 10 -60 is much larger than the trun- 
cation and round-off errors of the CNS results gained by means of the high-order 
Taylor expansion and data in 300-digit precision with a reasonable time step At, as 
illustrated above. So, the CNS provides us a convenient tool to study the transfer 
and evaluation of such kind of inherent micro-level uncertainty of initial conditions. 

The key point is that such an inherent micro-level uncertainty in the initial con- 
ditions finally leads to the huge, observable difference of orbits of the three bodies: 
it disrupts not only the elegant symmetry of the orbits but also the three-body sys- 
tem itself. Note that, Body 2 escapes with Body 3 in the case of 5 = +10~ 60 , but 
with Body 1 in the case of 5 = — 10~ 60 , to become a binary-body system! In nature, 
such kind of inherent micro-level uncertainty exists for each body at any time t > 0. 
So, from the physical viewpoint, the orbits of each body at large enough time are 
inherently unknown, i.e. random: given the same initial condition (in the viewpoint 
of physics), the orbits of the three-body system under consideration might be com- 
pletely different. For example, the three-body system might either have the elegant 
symmetry, or disrupt as the different binary-body systems, as shown in Figs. [9] to 11 



and Figs. [Eq to 14 , respectively 



Note that, from mathematical viewpoint, we can accurately simulate the orbits 
of three bodies in the interval < t < 1000. However, due to the inherent posi- 
tion uncertainty and the SDIC of chaos, orbits of the three bodies are random when 
t > 810, since the inherent position uncertainty transfers into macroscopic random- 
ness. Thus, from physical viewpoint, there exists the maximum predictable time T^ ax1 
beyond which the orbits of the three bodies are inherently random and thus can not 
be predictable in essence. Note that Z^ a3 , is determined by the inherent position un- 
certainty of the three bodies, which has nothing to do with the ability of human being. 
Therefore, long term "prediction" of chaotic dynamic system of the three bodies is 
mathematically possible, but has no physical meanings! 

Finally, to confirm our above conclusions, we further consider such a special case 
with the micro-level uncertainty of the initial position dv\ = 10 _60 (1, 1, 1), i.e. 

ri = (o,o,-i) + io- 60 (i,i,i). 

It is so tiny that, from the physical viewpoint mentioned above, the initial positions 
can be regarded as the same as those of the above-mentioned three cases. However, 
the corresponding orbits of the three bodies (in the time interval < t < 1200) 
obtained by means of the CNS with At = 10~ 3 , iV = 300 and the 30th-order Taylor 
expansion (M = 30) are quite different from those of the three cases: Body 2 first 
oscillates along a line on z = but departs from the regular orbit for large t to move 
along a complicated three-dimensional orbits, while Body 1 and Body 3 first move 



with the symmetry but lose it for large t, as shown in Figs. 15 to 17. However, it 
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is not clear whether any one of them might escape or not, i.e. the fate of the three- 
body system is unknown. Since such kind of micro-level uncertainty of position is 
inherent and unknown, given the same (from the physical viewpoint) initial positions 
of the three bodies, the orbits of the three-body system at large enough time is 
completely unknown. So, it has no physical meanings to talk about the accurate 
long-term prediction of orbits of the three-body system, because the orbits at large 
time (such as t > 1000) is inherently unknown/random and thus should be described 
by probability. This is quite similar to the motion of electron in an atom. It should 
be emphasized that, such kind of transfer from the inherent micro- level uncertainty 
to macroscopic randomness is essentially due to the SDIC of chaos, but has nothing 
to do with the ability of human being and Heisenberg uncertainty principle [4] . 

All of these reliable CNS results indicate that, due to the SDIC of chaos, such kind 
of inherent micro-level uncertainty of a star /planet might transfer into macroscopic 
randomness. This provides us an explanation for the macroscopic randomness of the 
universe, say, the inherent micro-level uncertainty might be an origin of the micro- 
scopic randomness, although it might be not the unique one. This might enrich and 
deepen our understandings about not only the three-body problem but also the chaos. 



4 Conclusions 

The famous three-body problem is investigated by means of a numerical approach 
with negligible numerical noises in a long enough time interval, namely the Clean Nu- 
merical Simulation (CNS). From physical viewpoints, position of any bodies contains 
inherent micro-level uncertainty. The evaluations of such kind of inherent micro-level 
uncertainty are accurately simulated by means of the CNS. Our reliable, very accu- 
rate CNS results indicate that the inherent micro-level uncertainty of position of a 
star/planet might transfer into macroscopic randomness. Thus, the inherent micro- 
level uncertainty of a body might be an origin of macroscopic randomness of the 
universe. In addition, from physical viewpoints, orbits of some three-body systems 
at large time are inherently random, and thus it has no physical meanings to talk 
about the accurate long-term prediction of the chaotic orbits. Note that such kind of 
uncertainty and randomness has nothing to do with the ability of human being and 
Heisenberg uncertainty principle 131. 

In this article, we introduce a new concept: the maximum predictable time T!^ az 
of chaotic dynamic systems with physical meanings. For the chaotic motions of the 
three body problem considered in this article, there exists the so-called maximum 
predictable time T^ ax m 810, beyond which the motion of three bodies is sensitive to 
the micro-level inherent uncertainty of position and thus becomes inherently random in 
physics. The so-called maximum predictable time of chaotic three bodies is determined 
by the inherent uncertainty of position and has nothing to do with the ability of human 
being. Note that, from the mathematically viewpoint, we can accurately calculate 
the orbits by means of the CNS in the interval < t < 1000. Thus, considering the 
physical inherent uncertainty of position, long term "prediction" of chaotic motion of 
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Figure 15: Orbit of Body 1 (0 < t < 1200) when dr ± = 10" 60 (1, 1, 1). 




Figure 16: Orbit of Body 2 (0 < t < 1200) when dr 1 = 10" 60 (1, 1, 1) 




Figure 17: Orbit of Body 3 (0 < t < 1200) when dr 1 = 10" 60 (1, 1, 1) 



the three bodies is mathematically possible but has no physical meanings. 

In summary, our rather accurate computations based on the CNS about the famous 
three-body problem illustrate that, the inherent micro-level uncertainty of positions 
of starts/planets might be one origin of macroscopic randomness of the universe. This 
might enrich our knowledge and deepen our understandings about not only the three- 
body problem but also chaos. Indeed, the reliable computations based on the CNS 
are helpful for us to understand the world better. 

Note that the computation ability of human being plays an important role in the 
development of chaotic dynamic systems. The finding of SDIC of chaos by Lorenz 
in 1963 is impossible without digit computer, although data used by Lorenz in his 
pioneering work is only in accuracy of 16-digits precision. So, the CNS with negligible 
numerical noises provides us a useful tool to understand chaos better. 



Finally, as reported by Sussman and Jack 14 15 , the motion of Pluto and even 
the solar system is chaotic with a time scale in the range of 3 to 30 million years. 
Thus, due to the SDIC of chaos and the micro-level inherent uncertainty of positions 
of planets, the solar system is in essence random. Note that such kind of randomness 
has nothing to do with the ability of human: the history of human being is indeed 
too short, compared to the time scale of such kind of macroscopic randomness. The 
determinism is only a concept of human being: considering the much shorter time- 
scale of the human being, one can still regard the solar system to be deterministic, 
even although it is random in essence. 
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